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We present first results from a full QCD calculation on the CP-PACS, comparing various actions at a -1 ~ lGeV 
and m^/m p w 0.7-0.9. We use the plaquette and a renormalization group improved action for the gluons, and 
the Wilson and the SW-Clover action for quarks. We find that significant improvements in the hadron spectrum 
results from improving the quarks, while the gluon improvement is required for a rotationally invariant static 
potential. An ongoing effort towards exploring the chiral limit in full QCD is described. 



1. Introduction 

With the progress in recent years of quenched 
simulations of QCD, deviations of the quenched 
hadron spectrum from experiment are becoming 
apparent. Most recently the precision quenched 
results from the CP-PACS show a systematic de- 
viation not only in the meson sector, but also in 
the spectrum of baryons [0 . Clearly the time has 
come to bolster efforts towards simulations of full 
QCD. 

Full QCD simulations are, however, extremely 
computer time consuming compared to those of 
quenched QCD. This leads us to consider im- 
proved actions for our simulation of full QCD. 

Studies of improved actions have been widely 
pursued in the last few years. However, a system- 
atic investigation of how various terms added to 
the gauge and quark actions, taken separately, af- 
fect light hadron observables has not been made 
in full QCD. In this article we report results of 
such a study carried out on the CP-PACS. 



* Talks presented by R. Burkhalter and T. Kaneko at 
Lattice97, Edinburgh, Scotland, 22-26 July 1997. 



2. Actions and simulation parameters 

We are carrying out a comparative study of 
the light hadron spectrum and static quark po- 
tential for four action combinations. For gluons 
we choose either the standard plaquette action 
(P) or an action (R) which was obtained by a 
renormalization group treatment Q and includes 
a rectangular Wilson loop. For the quarks we use 
either the Wilson action (W) or the Clover action 
(C) H). For the clover coefficient csw we com- 
pare: 

(a) the tree value csw = 1< 

(b) the meanfield (MF) improved value [§J csw = 
_P~ 3 / 4 with P the self-consistently determined 
plaquette average, and hence a different value of 
csw for each K. 

(c) a perturbative meanfield (pMF) improved 
value c SW = (1 - 0.8412 • /J" 1 )^ 3 / 4 , employed for 
the R action, where the one-loop value is used for 
the plaquette. 

We expect the extent of improvement to be 
clearer at a coarser lattice spacing. We there- 
fore attempt to tune the coupling constant (3 so 
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Table 1 

Simulation parameters. All the runs are on a 
12 3 x32 lattice, except (*), which is on a 16 3 x32 
lattice. For runs marked with (f ) the static quark 
potential is measured. Errors in the mass ratios 
are statistical but for a -1 they are systematic, 
taking into account uncertainties from the chiral 
extrapolation. 
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that the lattice spacing equals a -1 ~ lGeV. 

Our simulations are carried out for two flavors 
of quarks, mostly on a 12 3 x 32 lattice. We em- 
ploy the hybrid Monte Carlo algorithm to gen- 
erate full QCD configurations, generally at three 
values of K corresponding to m^/ra p rs 0.7-0.9. 
Details of the simulation parameters are given in 
Table |. For the inversion of the fcrmion matrix 
we first used the MR algorithm but switched later 
to BiCGStab. The molecular dynamics step size 
is chosen to yield an acceptance of 70-90%. 

We measure hadron propagators every 5 tra- 
jectories, using smeared sources and point sinks 
following the method of our quenched study 0]. 
Errors of masses are determined with a jackknifc 
analysis using a bin size of 5. The lattice spacing 
is fixed by the p meson mass extrapolated to the 
chiral limit through the form m p = m p ' + C2TO 2 . 
or +C2m^ +C3TO 3 . The errors for a" 1 quoted 
in Table [l] include the systematic errors in the ex- 
trapolation, estimated from the two fitting func- 
tions above, and by varying the fitting ranges. 

Measurements of the static quark potential are 
performed on about 100 configurations at those 
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PC 5.0 
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parameters marked with a t in Table [l| In these 
measurements we employ the smearing procedure 
proposed in Ref. || and extract the ground state 
potential in the same way as in Ref. || . 

3. Light hadron spectrum 

Our main results for the effect of improved ac- 
tions on hadron masses are displayed in Fig. |l|, 
in which the ratio rri]\r/m p is plotted as a func- 
tion of (tott/top) 2 for the four action combina- 
tions. The solid curve represents the well-known 
phenomenological mass formula . 

For the standard action combination PW (up- 
ward triangles), the ratios are well above the 
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Figure 1. mpf/m p as a function of (m T /m p ) 2 for 
four combinations of the action. Clover results 
are for the choice csw = MF. 
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Figure 2. Scaling behaviour of ra^jmp at fixed 
m^/mp as a function of m p a for various combi- 
nations of the action. The star on the left is the 
prediction of the phenomenological formula. 



phcnomenological curve. When we improve the 
gauge action (the RW case, shown by down- 
ward triangles), the data points come closer to 
the curve. By far the most conspicuous change, 
however, is observed when we introduce the clover 
term to the quark action. For both the PC (cir- 
cles) and RC (squares) cases, the data points lie 
much nearer to, or on top of, the phcnomcnolog- 
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Figure 3. m p and m« as a function of m 2 for 
different choices of csw- Dashed lines are chiral 
fits, cubic in m^, to data with csw = MF. 

ical curve. 

Looking only at Fig. [[], however, is mislead- 
ing because the runs for different action combi- 
nations are performed at slightly different lat- 
tice spacings. In Fig. || we therefore show the 
a-dependence. We interpolate the ratios m^lm p 
to a common, fixed ratio /m p = 0.8 and plot 
m]v/m p vs. m p in lattice units, with m p also 
interpolated to m w /m p = 0.8. Now it becomes 
clear that the better results for the RW case are 
only due to a finer lattice spacing and that the 
PW and RW results seem to lie almost on the 
same curve. The results for the PC and RC cases 
with the MF or pMF choice of csw also lie on 
a common curve, which, however, is much flat- 
ter. Somewhat less, but still quite substantial, 
improvement is found for the tree level csw- The 
same trend is seen also at m w /m p = 0.7 and for 
the ratio m&/m p ||. 

Changing the gluon and quark actions one at 
a time, we have been able to see clearly that the 
clover term plays a decisive role in improving the 
spectrum. In this regard, improving the gluon 
action has much less effect. 

Another interesting feature in our hadron mass 
data is that they exhibit a negative curvature in 
terms of m 2 towards the chiral limit. In Fig. || 
we show this effect for the p meson and nucleon 
masses using the RW and RC action combina- 
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Figure 4. "Speed of light" as a function of the 
meson mass for four action combinations. Clover 
results are for the choice csw — MF 

tions and several choices of csw- We also note 
that for the Clover action the p meson masses are 
shifted to a higher value compared to those for 
the Wilson action. For the nucleon, however, we 
sec no difference between the Wilson and Clover 
cases. This confirms an earlier observation using 
staggered dynamical fermions|^]. 

4. Dispersion relation 

Improving the action is expected to lead to a 
dispersion relation closer to that in the contin- 
uum. As a measure for the closeness to the con- 
tinuum we define the "speed of light" c through 
c 2 = dE{p) 2 /dp 2 \ p=0 . Values of c are extracted by 
a linear fit of E(p) 2 to p 2 , using the three lowest 
momenta. In the continuum c equals one for all 
values of the rest mass. On the lattice, however, 
c is only expected to go to one for vanishing rest 
mass. For nonzero rest mass it is generally differ- 
ent from one. An improved action is expected to 
have a "speed of light" closer to one for nonzero 
rest masses. 

Our results for c are plotted in Fig. |4| An im- 
provement is most clearly seen in a comparison 
of the PW and RC action combinations. We 
are not able to tell whether the improvement is 
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Figure 5. Static quark potential for four action 
combinations at m^/nip ~ 0.8. Scales are nor- 
malized by the lattice spacing determined from 
TO p in the chiral limit. Different symbols cor- 
respond to potential data measured in differ- 
ent spatial directions along the vector indicated 
in the figure. Clover results are for the choice 
c sw = MF 

due more to the gauge or quark action, however, 
since the RW and PC cases both also show signs 
of similar improved behavior within the present 
statistics. 

5. Static quark potential 

5.1. Rotational symmetry 

The static quark potential provides a good 
measure of the effectiveness of improving the ac- 
tion in restoring the rotational symmetry. In the 
quenched case it is known that the symmetry is 
well restored by improving the gauge action ^|[l0[. 
Our typical results for full QCD taken for a quark 
mass corresponding to m^/nip w 0.8 are plotted 
in Fig. for each action combination. 

We see quite clearly in Fig. || that the symme- 
try is remarkably well restored by improving only 
the gauge action also in full QCD. The effect of 
quark action improvement on the restoration of 
the symmetry is not clear in our data. 

This may be explained by the fact that the 
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Figure 6. Lattice spacing a determined by 
the string tension a and the p meson mass as 
a function of (m^a) 2 for the four action combi- 
nations. Clover data are for the choice csw — 
MF. Filled circles are obtained by identifying 
era 2 = (440McV) 2 , and open circles from m p a = 
770MeV at each values of K and in the chiral 
limit. 

clover improvement does not remove terms of 
the Wilson action breaking rotational symme- 
try, while the R action for gluons has much 
reduced coefficients for the rotationally variant 
terms compared to the plaquette action. 

5.2. Matching of lattice scale 

Another interesting question one can address 
with data for the static potential is whether the 
lattice scale determined from the string tension a 
is consistent with that from the p meson mass. 

We extract the string tension assuming the 
form V(r) = Vq — a/r + a-r. The Coulomb co- 
efficient a is ill-determined with our data, which 
starts at a coarse distance of r min = a « 0.2 fm. 
We therefore fix a = 0.3, where x 2 /dof is gen- 
erally the smallest, and use the shift of a over 
the range a = 0.2 — 0.4 as a systematic error. 
This fitting is repeated over several fitting ranges, 
and the average of the fitted results is taken as 
the central value of <r, while the variance is in- 
cluded into the systematic error. For estimating 
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Figure 7. Static quark potential obtained with 
RC action with csw = pMF on an 16 3 x32 lat- 
tice. 

the lattice spacing we employ the phenomenolog- 
ical value a = (440MeV) 2 . In Fig. § we plot the 
lattice spacing obtained from the string tension 
and the p meson mass for each value of K for all 
four action combinations. For the p meson mass 
the value in the chiral limit is also shown. 

Let us emphasize that the consistency of the 
two determinations of the lattice scale is expected 
only in the chiral limit and sufficiently close to the 
continuum limit. From this point of view, a sig- 
nificant discrepancy in the chiral limit between 
the two determinations, as observed in the PW 
and RW results, shows the presence of large scal- 
ing violating effects, originating from the Wilson 
quark action, at a -1 w 1 GcV. We expect the 
discrepancy to disappear closer to the continuum 
limit, and this is supported by the results ob- 
tained at $ = 5.5 with a' 1 w 2 GeV of Ref. 0. 

In contrast, the two estimates of the lattice 
spacing converge well to a consistent result for 
the PC and RC cases. This shows that the clover 
term helps to improve the consistency of the two 
determinations of the lattice spacing already at 
a- 1 w l.OGeV. 

5.3. Sea quark effects in the potential 

Our results for the potential, plotted in Fig. ||, 
and covering distances up to r w 1 fm, are quite 
linear in r and do not show any sign of string 
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Figure 8. Results on a 16 3 x 32 lattice for 
RC action at (3 = 1.9. Diamonds are exper- 
imental points corresponding to N(940)//?(770), 
A(1232)/p(770) and ft(1672)/0(lO2O). 

breaking via qq pair creation. In order to search 
for this effect, we employ one of the RC runs 
carried out on a 16 3 x 32 lattice at = 1.9 with 
which we can measure the potential up to r ~ 
2fm. The results are shown in Fig. |?], and still 
do not show any sign of deviation from a linear 
behaviour in r. 

We find, however, that the maximum overlap 
of the ground state with the smeared Wilson loop 
operator decreases as r increases and falls down to 
about 50% at r ~ 2.0 fm for our full QCD case, 
while effectively a 100% overlap is obtained for 
the quenched case, independent of the distance 
r. This may be an indication that the ground 
state develops towards larger r a mixing, absent 
in quenched QCD, with states containing qq pairs. 

6. Towards the chiral limit 

A very interesting and important issue in simu- 
lations in full QCD is how close one can approach 
the chiral limit. Of the four action combinations 
we have examined, the RC combination exhibits 
the most improved behaviour with regard to the 
scaling of mass ratios and the rotational symme- 
try of the potential. We therefore choose this ac- 
tion combination for a feasibility study towards 
the chiral limit. The runs are made at j3 = 1.9 on 
a 16 3 x32 lattice with a spatial size of about 2.6 fm, 



and a quark mass corresponding to m^/m p as 
small as ~ 0.4 is explored (see Table ||). Results 
for the mass ratios obtained with 80-260 configu- 
rations for each quark mass, which includes statis- 
tics accumulated after the conference (especially 
for the lightest quark), are shown in Fig. [| 

These results point towards the possiblity that, 
by improving both the fermion and the gauge 
action, continuum results in full QCD at small 
quark masses are achievable using the presently 
available computing power. 
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